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FOREWORD 

This report summarizes the work conducted by McDonnell Douglas Astronautics 
Company-St. Louis (MDAC-St. Louis) for the NASA Langley Research Center under 
Contract NAS1-12436. Mr. James C. Dunavant was the technical monitor for the 
study. Mr. H. J. Fivel (MDAC-St. Louis) made a significant contribution in 
developing the computer coding used in the assimilation and analysis of the wind 
tunnel data presented in this report. Dr. R. T. Krieger, Mr. J. K. Lehman, and 
Dr. R. M. Laurenson contributed to the aerodynamic, material, and structural 
dynamic aspects of the study. Mr. Darrell Weber assisted in formulating and 
coding the heating models used in the thermal flight assessment. 
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1 . 0 SUMMARY 


A study has been conducted to investigate the feasibility of using corrugated 
panels as a thermal protection system (TPS) for an advanced space transportation 
vehicle. The study consisted of two major tasks: (1) development of improved 
correlations for wind tunnel heat transfer and pressure data to yield design 
techniques, and (2) application of the design techniques to determine if corrugated 
panels have application on future aerospace vehicles. A si ngle-stage-to-orbi t 
(SSTO) vehicle was used in this study to assess advantages and aerothermodynamic 
penalties associated with use of such panels. 

In the correlation task, experimental turbulent heat transfer and pressure 
data obtained on corrugation roughened surfaces during wind tunnel testing were 
analyzed and compared with flat plate data. The correlations and data comparisons 
include the effects of a large range of geometric, inviscid flow, internal boundary 
layer, and bulk boundary layer parameters in supersonic and hypersonic flow. 

The formulated wind tunnel correlations were used in the second task to 
evaluate the increased surface temperatures and added insulation requirements 
associated with the use of corrugated panels on an advanced space transportation 
system. The thermal evaluation considered the effects of panel location, material 
selection, and lateral heat conduction due to the large structural temperature 
gradients normal to the corrugation axes. The second task also included an 
analysis of the impact of corrugated panels on reentry vehicle aerodynamics and 
an evaluation of the flutter characteristics associated with corrugated panels 
on a SSTO in the supersonic and hypersonic range. 


1-1 



2.0 INTRODUCTION 

The thermal protection system (TPS) is an important factor in the design of 
advanced space transportation systems. The space systems of current interest 
range from reusable, ballistic-entry systems to single-stage-to-orbit (SSTO), 
lifting-entry vehicles. Interest in the present study is focused on the SSTO 
which will be a winged reentry vehicle anticipated to be operational in the 1990's 
and beyond. The TPS for the SSTO must be reusable, rugged, and forgiving of 
abnormal operating conditions. The Reusable Surface Insulation (RSI) being used 
for the Space Shuttle TPS is fragile and may not be applicable for the SSTO due 
to stringent service life requirements. 

A reusable, non-fragile TPS which is being studied for advanced space 
transportaion systems is composed of radiation-cooled, metallic panels (Ref. 1). 

The panels are stiffened by corrugations or waves which protrude above the equivalent 
smooth surface of the vehicle and distort the boundary layer flow during flight. 

The overall effect of this flow distortion on the vehicle performance character- 
istics is a significant technical consideration which must be considered during 
the development of corrugated panel TPS. The emphasis of this study is focused on 
the effect of surface corrugations on the thermal characteristics of an advanced 
reentry vehicle. The effect of the corrugations on the flight performance is also 
considered, but in less detail. 

A meaningful investigation of the boundary layer distortion due to surface 
corrugations must address the situation for which the surface corrugations are 
deeply submerged in the thick turbulent boundary layers which will develop on 
large vehicles such as the SSTO. Several investigators have studied the influence 
of surface waves on turbulent boundary layers (Refs. 2-7). However, these studies 
did not include data for which the boundary layer was orders of magnitude thicker 
than the corrugation height. The first detailed investigation which considered 
heat transfer to corrugations deeply submerged in a thick turbulent boundary layer 
was reported by Brandon, Masek, and Dunavant (Refs. 8 and 9). . During that study, 
questions arose regarding the influence of several geometric parameters on the 
boundary layer distortion. As a result, NASA/LRC conducted additional wind tunnel 
experiments to study the effect of the wave shape on the boundary layer. 

In this report the additional data for the effect of wave shape are analyzed, 
correlated, and compared with the previous data obtained in thick and thin tur- 
bulent boundary layers. The correlations developed in Refs. 8 and 9 are improved 
by including data for additional wave configurations and flow conditions and by 


2-1 


ORIGINAL PAGE IS 
OF POOR QUALITY 



updating the boundary layer parameters to conform with recent wind tunnel experiments 
The results of the wind tunnel data analysis are then used to assess the feasibility 
of using corrugated panels on a SSTO vehicle. In the flight assessment, available 
engineering design procedures are used to compute the inviscid and boundary layer 
flow for a given flight trajectory. The flight vehicle flow field calculations are 
then used in conjunction with the wind tunnel correlations to examine the heating 
penalties associated with corrugated panels at several locations on the vehicle. 

The wind tunnel data and flight flow field calculations are also used to deduce the 
reduction in aerodynamic performance and structural integrity due to the corrugations 
It should be noted that an engineering solution was chosen for the flight flow 
field calculations because it was felt that a more exact solution was not warranted. 
It was thought that if the results of an engineering solution showed that 
corrugated panels appear feasible, a more exact analysis could then be justified 
in a benchmark assessment before conducting a flight test program. 
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3.0 MODEL AND INSTRUMENTATION 

Three corrugation roughened flat panels were tested in this program. The 
panels were similar to those fabricated and tested in Ref. 8 . All three panels 
were 50.8 cm x 50.8 cm x 2.54 cm and were fabricated from nominally 0.051 cm thick 
347 stainless steel. As shown in Figure 3-1, the panels were designated -1, -2, 
and -3 and had peak amplitudes of 1.06, .525, and .405 cm, respecti vely. The 
wavelength was 7.32 cm for panels -1 and -2, and 3.66 for panel -3. Cross-section 
shapes for panels -1, -2, and -3 are compared with the shapes previously tested 
(designated -4 and -5) in Ref. 8 and with the configuration tested by Sawyer 
in Ref. 10 (designated panel -6).* All the shapes shown in this figure are for 
sections normal to the corrugation axes. Further details of the overall panel 
shape for -1, -2, and -3 are presented in Figures 3-2 and 3-3. For all three 
panels the amplitude was constant over a 25.4 cm section in the center, and tapered 
to zero near the panel edge. Panels -1 and -2 contained 6 corrugations (wave cycles) 
and panel -3 contained 12 corrugations which ran the entire length of the model. 

The corrugation sections were constructed of circular arcs connected by straight 
line segments, as shown in Figure 3-4. Also shown in this figure is the nomen- 
clature associated with the wave construction. 

The effective wave shapes for panels -1, -2, -3, and -6 were computed 
as a function of the flow angle relative to the corrugations, 0, using elliptical 
integrals. It was apparent that the wave shapes were greatly influenced by the 
flow angle. For example, the maximum surface deflection angle varied from 15.1° 
for 0 = 15° to 43.9° for 0 = 90° for panel -1, from 4.6° for 0 = 15° to 17.2° 
for 0 = 90° for panel -2, and from 7.4° for 0 = 15° to 26.6° for 0 = 90° for 
panel -3. The maximum surface deflection angle for Panel -6 (Sawyer panel) 
varied from 6.9° for 0 = 15° to 30.1° for 0 = 90°. 

Panels -1, -2, and -3 were instrumented with thermocouples and pressure 
orifices. These panels were instrumented in much more detail than those tested 
in Ref. 8 so that a much better definition of the heating and pressure distributions 
was obtained. With the waves normal to the flow direction (configuration termed 
0 = 90°), corrugations in the center of the panel were instrumented in the flow 


*To facilitate reference to the corrugation shapes as they are discussed in the text, 
Figure 3-1 is reproduced as a foldout page in Appendix A. 

**Note 0 = 90° indicates the flow is normal to the corrugation axis. 
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FIGURE 3-1 CORRUGATION SHAPES INVESTIGATED 
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NOMINAL VALUES GIVEN ARE NONDIMENSIONALIZED BY 2.54 cm 
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FIGURE 3-4 BASIC WAVE CONSTRUCTION AND NOMENCLATURE 
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direction. The first complete cycle and the last complete wave cycle were also 
instrumented. However, only the data in the center of the panel are of interest 
to the present study since they have previously been shown to be representative 
of all the waves in very thick turbulent boundary layers (neglecting edge effects). 

Details of the geometry and instrumentation for panels -4 and -5, and -6 
are discussed in Ref. 8 and 10, respectively. 
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4.0 TEST PROGRAM 


Panels -1, -2, and -3 were tested in the turbulent wall boundary layer 
of the Langley Continuous Flow Hypersonic Wind Tunnel (CFHT) which has a 
freestream Mach number of 10.3. The wall boundary layer provided simulation 
of the thick boundary layers occurring on a flight vehicle. The procedure 
used in the test program was the same as that described in Ref. 8. The panels 
were attached to an adapter plate which was mounted flush with the wind tunnel 
wall, as shown schematically in Figure 4-1. A schematic of the corrugation 
angle, is also shown in this figure. For the CFHT tests, the corrugation 
orientation angle relative to the freestream direction could be varied from 
0° to 15° and from 75° to 90°. Angles between 15° and 75° could not be ob- 
tained due to the arrangement of the tunnel injection mechanism. Panels -1, 

-2, and -3 were not tested in the Unitary Plan Wind Tunnel (UPWT) during the 
present test program. 

Heating and pressure data obtained on panels -4 and -5 and previously ! 
reported in Refs. 8 and 9 are also considered in this study. These two panels 
were not further tested for the present study. In this study these previously 
obtained data are used in conjunction with improved boundary layer calculations. 
The data obtained by Sawyer (Ref. 10) are also used in this study. The Sawyer 
data were obtained on a corrugated panel configuration which was tested in the 
UPWT at freestream Mach numbers of 2.5 and 4.5. 

The flow and panel geometric conditions for the data considered in the 
present study are listed in Table 4-1. Flat plate data were also obtained for 
all the Mach number and Reynolds number combinations. The boundary layer para- 
meters for the turbulent wall boundary layers are listed in Table 4- II. These 
values differ from those published in Ref . 8 which were computed assuming that 
the freestream static pressure was constant through the boundary layer. Sub- 
sequent wall pressure measurements in the CFHT showed that this was not the 
case. Hence, all the CFHT calculations were repeated using measured pitot and 
total temperature profiles and an assumed linear static pressure distribution 
from the freestream static pressure to the measured wall pressure. The com- 
puted boundary layer profiles for the CFHT are shown in Figures 4-2 to 4-4 
as a function of the freestream Reynolds number. The boundary layer data 
measured in Refs. 11 and 12 were used to obtain the UPWT boundary layer para- 
meters. For both the CFHT and UPWT, the laminar sublayer thickness was com- 
puted from Beckwith's correlation (Ref, 13). 
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FIGURE 4-3 STATIC TEMPERATURE PROFILE IN CFHT WALL BOUNDARY LAYER 
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5.0 EXPERIMENTAL HEATING AND PRESSURE DISTRIBUTIONS 


The heating and pressure distributions measured on the present panels (-1,-2, 
and -3) are used in conjunction with the previously obtained data to deduce the 
effect of the flow angle relative to the corrugations (corrugation angle), Mach 
number, Reynolds number, wave length, wave amplitude, wave surface radius, and wave 
scaling. Previous tests in thick turbulent boundary layers have shown that the 
heating and pressure distributions are relatively insensitive to the location of 
the wave in the train. Hence, the present study is concerned only with data in 
the center of the panels. Distributions over a wave in the center are representative 
of those over all the waves (neglecting edge effects). All the heating and pressure 
values appearing in this report are nondimensional ized with respect to local measured 
flat plate values. This is necessary due to spanwise gradients which existed on 
the wall of the wind tunnel. Analysis has shown that the effects of the nonuniform 
boundary layer are eliminated by referencing the corrugated data to the local flat 
plate values. 

Flow Angle (Corrugation Angle) Effect - Heat transfer and pressure distributions 

r 

for hypersonic flow (M^ = 10.3, Re M /M = 6.6 X 10 ) are presented in Figures 5-1 to 
5-4 for panels -1 and -3 at flow angles of 0°, 15°, 75°, and 90°. Distributions 
for supersonic flow (M = 2.4, Re^/M = 10 X 10 ) for panel -5 at flow angles of 

0°, 15°, 30°, 60°, and 90° are presented in Figures 5-5 and 5-6. For all flow 

angles heating is much more affected by the presence of corrugations than is 

pressure. The largest heating and pressure increases occur when the flow is normal 

to the corrugations. These trends are similar to those observed in Ref. 8. Oil flow 
visualization indicated that these distributions were caused by flow separation in 
the valley and subsequent reattachment on the following wave. The greatest increase 
in heating over smooth wall values occurs as the flow angle increases from 0° to 
15°. However, the relative changes are more pronounced for panels -1 and -3. The 
increase in heating occurring from 0 = 75° to 90° is quite small for all three panels. 
The greatest increase in pressure over smooth wall values occurs over the range in 
flow angle from 15° to 75°. The increase in pressure from 0 = 75° to 90° is quite 
small for panels -1 and -3. These heating and pressure distributions could be 
affected by flow angularity in the wall boundary layer which could be a function of 
the distance normal to the wall. 

Mach Number Effect - Changes in the heating and pressure distributions on a 
corrugated panel due to increases in the local free-stream Mach number (while the 
free-stream Reynolds number is held constant) can be seen by comparing Figures 5-5 
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FIGURE 5-1 EFFECT OF FLOW ANGLE ON HEATING DISTRIBUTION, M = 10.3 
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FIGURE 5-2 EFFECT OF FLOW ANGLE ON PRESSURE DISTRIBUTION, M =10.3 
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FIGURE 5-4 EFFECT OF FLOW ANGLE ON 





FIGURE 5-5 EFFECT OF FLOW ANGLE ON HEATING DISTRIBUTION, PANEL -6 
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FIGURE 5-6 EFFECT OF FLOW ANGLE ON PRESSURE DISTRIBUTION 5 PANEL -6 



and 5-6 with Figures 5-7 and 5-8 in which data are presented for flow angles of 
0°, 15°, 30°, 60°, and 90° for M =2.4 and 4.5 for a constant free-stream 
Reynolds number per meter of 10 X 10°. The Mach number has a larger effect on 
the heating than on pressure with the importance being more pronounced at the 
higher flow angles. For both pressure and heating the increases over smooth 
wall values become greater with increasing Mach number. 

Reynolds Number Effect - Previous experiments on the influence of Reynolds number 
on corrugated panels in thick turbulent boundary layers in hypersonic flow were 
conducted at R 6co /M = 1.3 X 10®, 3.3 X 10®, and 6.6 X 10®. Results showed that peak 
heating and pressure increased with Reynolds number. Similar results were obtained 
in this study for the same Reynolds number range, as shown in Figures 5-9 and 5-10 
which are for free-stream flow normal to the waves of panel -1 at a free-stream 
Mach number of 10.3 and Re - 1.3, 3.3, 6.6, and 8.2 X 10®. Increasing R P in- 
creases the heating and pressure except for changes from 6.6 X 10® to 8.2 X 10® 
over which range the peak heating slightly decreases. Identical results were 
obtained for panels -2 and -3 except the heating slightly increased from Re^/M = 

6.6 to 8.2 X 10®. The trends may be somewhat influenced by flow angularity in 
the boundary layer which could be a function of both Reynolds number and distance 
from the wall . 

For supersonic flow (M^ = 2.4) over panel -6, the effect of Reynolds number 
can be seen by comparing the distributions in Figures 5-11 and 5-12 with Figures 

5-5 and 5-6 for flow angles of 15°, 30°, 60°, and 90°. As the Reynolds number 

6 6 o 

increases from 3.3 X 10 to 10 XI 0, the heating decreases except for 0=0 and 

15° for which the heating slightly increases. The peak pressure decreases slightly 

with increasing Reynolds number for 0 = 30°, 60°, and 90° at approximately the same 

X/L location on the wave. For 0 = 15°, the peak pressure also decreases slightly 

with increasing Reynolds number. However, the peak pressure location for this low 

flow angle has moved downstream near the crest of the wave for the lower Reynolds 

number case. 

Wave Length Effect - The wave length is one of the primary parameters considered 
in the design of corrugated panels. Until the present study, the effect of the 
wave length on the heating distributions in thick turbulent boundary layers had 
not been investigated. In order to investigate the wave length effect, panel -2 
was constructed and tested and the results are compared with those for panel -4 
which has approximately the same wave height as panel -2 and 1/2 the wave length. 
The heating was found to be relatively insensitive to the wave length for all the 
flow conditions considered in this study, as demonstrated in Figures 5-13 and 5-14. 
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FIGURE 5-8 EFFECT OF FLOW ANGLE ON PRESSURE DISTRIBUTION, PANEL -6 
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FIGURE 5-10 EFFECT OF REYNOLDS NUMBER ON PRESSURE DISTRIBUTION, M = 10.3 
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FIGURE 5-11 EFFECT OF FLOW ANGLE ON HEATING DISTRIBUTION, PANEL -6 







FIGURE 5-12 EFFECT OF FLOW ANGLE ON PRESSURE DISTRIBUTION, PANEL -6 
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FIGURE 5-13 EFFECT OF WAVE LENGTH ON HEATING DISTRIBUTION, M = 10.3 
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FIGURE 5-14 EFFECT OF WAVE LENGTH ON HEATING DISTRIBUTION, M =10.3 
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The heating distributions on both panels are almost identical for the three Reynolds 
numbers shown. For the higher Reynolds number (R eoo /M = 6.6 X 10 6 ), the average 
heating for the larger wave length panel is somewhat greater due to the increased 
heating along the compression surface of the wave. The present results on the effect 
of wave length can be very important in the design of corrugated heat shields for 
application in very thick boundary layers since the wave length can be significantly 
changed without a corresponding heating penalty. 

Wave Amplitude Effect - Previous tests on the effect of corrugation height (wave 
amplitude) in thick turbulent boundary layers yielded results which were radically 
different than those for thin boundary layers. Heat transfer and pressure were 
found to be relatively insensitive to the wave amplitude for all flow conditions 
investigated. As a consequence of the previous results, panel -1 was constructed 
and tested to bridge the gap in the data between the very thick and thin boundary 
layer data. Panel -1 has a wave height which is almost twice that of panel -4, 
and a wave length which is twice that of panel -4. The heating on panels -1 and 
-4 are compared in Figures 5-15 and 5-16 for R Gco /M = 1.3 and 6.6 X 10®. Unlike the 
previous results in thick boundary layers, the present results show that the wave 
height can significantly affect the peak heating on corrugated panels if the 
corrugation protrudes far enough out into the boundary layer. Peak heating data 
for panel -5 which had an amplitude about 1/2 that of -4 was very similar to data 
for panel -4. However, peak heating data for panel -4 which has an amplitude about 
1/2 that of -1 is well below the data for panel -1. Comparison of the data for 
panels -2 and -4 shows that the wave length has a small influence on the peak 
heating, and hence the large difference between the data on panels -1 and -4 is 
attributed to the change in the wave height. The curves faired through the data 
for panel -4 ane slightly different from those presented in Ref. 8. The present 
curves were generated by a computer plot routine which used a spline fit, while 
the previous distributions were hand faired. 
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FIGURE 5-15 EFFECT OF WAVE HEIGHT ON HEATING DISTRIBUTION, M = 10.3 
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FIGURE 5-16 EFFECT OF WAVE HEIGHT ON HEATING DISTRIBUTION, M = 10.3 
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Wave Radius Effect - The panels of the present study have shapes constructed of 
circular arcs connected by straight lines as shown by the sketch in Figure 3-4. 

For a given wavelength and wave height, there is some leeway in choosing the arc 
radius in the panel design. Experiments were conducted with panel -3 and the 
results are compared with data for panel -4 to determine the importance of changing 
the surface radius by a factor of 2. Panel -3 has a surface radius equal to .457cm 
which is 1/2 the radius of panel -4. Both panels have the same wave length (3.66cm). 
Panel -3 was designed to have the same wave height as panel -4. However, due to 
spring back in the metal during construction, panel -3 has a wave height of .406cm 
compared to .61cm for panel -4. Heating on the two panels is quite similar con- 
sidering the surface radii are different by a factor of two, as shown in Figures 
5-17 and 5-18. The peak and average heating are slightly higher on panel -3 which 
has the smaller wave radius. Similar results were found for other flow angles at 
Mach 10.3. The increased heating on the sharper wave probably would have been more 
pronounced if its wave height had been exactly equal to that of panel -4. 

Wave Scaling - The shape of panel -1 is scaled up from panel -4 by a factor of 
about 2. (Again, the wave height is not quite to scale due to the spring back 
problem encountered during construction.) The two wave shapes are geometrically 
similar in that the wave shape for panel -1 can be obtained by uniformly increas- 
ing all the dimensions of panel -4, and the two profiles have identical surface 
deflection angles along the waves. However, comparison of the data obtained on 
panels -1 and -4 (Figures 5-15 and 5-16) shows that the heating distributions are 
not similar in that the distribution on panel -4 cannot be uniformly increased to 
obtain the distribution on panel -1. Therefore, a known heating distribution on 
a given panel shape cannot be used to determine the heating distribution on a 
geometrically similar panel by uniformly changing the heating on the compression 
and expansion surfaces by a constant value. The heating on the first panel can 
only serve as a guide to the distribution which can be expected on the .. ... 
second panel. Undoubtedly, a complex function of geometric, inviscid flow, 
internal boundary layer, and bulk boundary layer parameters relates the heating 
distributions on different panels whether or not the panels are geometrically 
similar. This type of function was not considered in this present study for the 
heating distributions. However, such a function was explored for maximum heating 
and is discussed in the correlation section of this report. 
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FIGURE 5-17 EFFECT OF WAVE RADIUS ON HEATING DISTRIBUTION, M = 10.3 
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FIGURE 5-18 EFFECT OF WAVE RADIUS ON HEATING DISTRIBUTION, M = 10.3 
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6.0 CORRELATION OF WIND TUNNEL DATA 

The maximum and average wave heating measured during the study on panels -1, 

-2, and -3 are summarized in Tables 6-1 and 6-1 1. The pressure drag coefficients 
measured on the panels are given in Table 6-111. The maximum heating values re- 
present the highest measured values. The average values were obtained by integrat- 
ing the heating distributions, and the pressure drag was computed by integrating 
the measured pressure distributions. The drag data for hypersonic flow are not 
correlated herein because the values were extremely low and the effect of the in- 
creased drag due to the corrugations would be insignificant in the hypersonic range. 
The maximum and average pressures measured on the panels are listed in Tables 6- IV 
and 6-V. Data obtained by Sawyer, Ref. 10, on panel -6 are also presented in Tables 
6-1 > 6- II, 6- IV, and 6-V and Figures 8-3 and 8-4 for comparison. 

For supersonic flow over panel -6, the average heating values are always 
greater than the corresponding flat plate value. For hypersonic flow, panels -2 
and -3 have average values less than flat plate values for some conditions. How- 
ever, the average heating for the largest amplitude panel (e = 1.06 cm) are all 
greater than the flat plate values. The changes in average heating with increasing 
corrugation angle are not consistent with the maximum heating trends in that average 
heating actually decreases with increasing flow angle for many conditions. 

Correlation of Maximum Pressure 

A detailed description of the flow field occurring over a wave of a corrugated 
panel is beyond the scope of the present investigation. Such a description would 
require a turbulent mixing solution with imbedded shocks coupled with a technique 
that predicts flow separation and the corresponding reattachment region. However, 
some effort was directed toward developing an understanding of the flow field over 
a wave using an elementary flow analysis. 

A simple flow model was investigated to determine if the maximum surface pres- 
sure could be predicted. In this model the local internal flow* at the wave height 
is assumed to expand isentropically over the expansion surface, pass through a 
separation shock (assumed to be a normal shock), and then be compressed isentropically 
at the stagnation point on the compression surface. The total pressure behind the 
imbedded shock would then represent the maximum surface pressure which could be 
recovered on the compression surface. Two cases were considered for which the 
separation and reattachment points were determined from oil flows. For both cases 
the flow was normal to the corrugations. The first case was supersonic flow 
(M = 2.4) over panel -6 and the second case was hypersonic flow (M ro = 10.3) over 

*M e is the Mach number at the wave maximum height determined in the flat plate 
boundary layer surveys. 
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TABLE 6-1 MAXIMUM WAVE HEATING 


Re /M 

00 

M 

oc 


<j> (DEG) 

PANEL -1 
(e = 1 . 06cm) 

h MAX /h fp 

PANEL -2 
(e = .525cm) 

h MAX /h fp 

PANEL -3 
(e = .405cm) 

h MAX /h fp 

I . 3x1 0 6 

10 

.3 

0 

1.45 

1.18 

1.32 





15 

3.66 

1.51 

1.57 

\ 

i 


i 

75 

3.33 

1.50 

1.88 

1 




90 

3.35 

1.70 

1.80 

3.3x1 0 6 



0 

1.36 

1.01 

1.12 





15 

3.44 

1.55 

1.69 





75 

4.36 

2.11 

2.09 


1 



90 

4.57 

2.15 

2.40 

6.6x10® 



0 

1.31 

.874 

1.00 





15 

2.63 

1.47 

1.92 





75 

4.54 

2.40 

2.61 

> 

f 



90 

4.63 

2.56 

2.90 

8.2> 

<10 6 



0 

1.31 

.870 

.981 





15 

2.43 

1.37 

1.90 



\ 

/ 

75 

4.38 

2.34 

2.70 

\ 

f 


f 

90 

4.41 

2.46 

2.98 


jeXshI 

3.3xl0 6 .79 2.4 


10. xl 0 6 


4.5 


y 


t 


(DEG) 

PANEL -6 
W h fp 

0 

1.26 

15 

1 .41 

30 

1.92 

60 

2.32 

90 

2.52 

0 

1.31 

15 

1 .48 

30 

1.84 

60 

2.10 

90 

2.30 

0 

1.46 

15 

1 .72 

30 

2.20 

60 

2.56 

90 

3.04 
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TABLE 6- I I AVERAGE WAVE HEATING 


PANEL -1 
(e = 1.06cm) 

h av /h fp 


PANEL -2 

(e = .525cm) 

h /h f 
av' fp 


h av /h fp 


1 . 3x1 0 6 10.3 



0 

1.06 

15 

1.85 

75 

1.20 

90 

1.25 

0 

1.16 

15 

1.81 

75 

1.52 

90 

1.62 

0 

1.10 

15 

1.53 

75 

1.68 

90 

1.72 

0 

1.10 

15 

1.42 

75 

1.69 

90 

1.68 


.920 

.890 

.741 

.836 

.839 

.912 

.893 

.922 

.792 

.903 

1.13 

1.22 

.798 

.862 

1.15 

1.26 


.966 

.954 

.917 

.912 

.936 

.925 

.891 

1.02 

.918 

1.04 

1.09 

1.26 


• 1 
1.03 
1.19 
1.34 


W M 

e(cm) 

3.3xl0 6 

.79 



(DEG) 

Panel -6 
h av /h fp 

0 

1.10 

15 

1.13 

30 

1.37 

60 

1.50 

90 

1 .62 

0 

1.16 

15 

1.21 

30 

1.37 

60 

1.53 

90 

1.61 

0 

1.24 

15 

1.20 

30 

1.43 

60 

1.44 

90 

1.69 
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TABLE 6-1 I I 


PRESSURE DRAG COEFFICIENTS 
TOR CFHT DATA 


panel! 

Re /M = 1.2 

00 

1 x 10 6 

3.3 x 10 6 

6.6 x 10 6 

8.2 x 10 6 

* = 15° 

3.57 x 

10" 5 

1.82 x 10" 5 

1.60 x 10~ 5 

1.55 X 10“ 5 

75° 

.00182 

• .00126 

.00120 

.00118 

90° 

.00128 

.00216 

.00143 . 

.00141 

PANEL 2 
* = 15° 

7.55 x 

10“ 7 

-7.66 x 10~ 7 

-1.10 x 10“ 6 

-1.15 x 10" 6 

O 

in 

2.08 x 

10’ 4 

1.62 x 10~ 4 

2.35 x 10" 4 

2.72 x 10" 4 

90° 

2.09 x 

10- 4 

2.63 x 10” 4 

3.58 x 10" 4 

3.94 x 10" 4 

PANEL 3 

<p = 15° 

4.40 x 

10' 6 

1.27 x 10" 6 

9.63 x 10“ 7 

8.20 x 10" 7 

75 

4.72 x 

10- 4 

3.18 x 10" 4 

3.68 x 10" 4 

4.19 x 10" 4 

0 

o 

<y> 

3.91 x 

10" 4 

4.73 x 10' 4 

5.41 x 10" 4 

5.77 x 10' 4 

c d = 

f pres. 

1 II 2 

2 P <A 

(A) 

I Based 

on wavelength in 

freestream flow 

di recti on 
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TABLE 6- IV MAXIMUM WAVE PRESSURE 


PANEL -1 PANEL -2 PANEL -3 


Re /M 

CO 

M 

00 


<p (DEG) 

(e = 1 .06cm) 
P max^ P fp 

(e = ,525cm) 

P /P f 
max fp 

(e = .405cm) 
P max^fp 

1.3xl0 6 

10.3 

0 

1.03 


1.19 

1.09 





15 

1.28 


1.21 

1.27 





75 

1.43 


1.42 

1.43 


f 



90 

1.19 


1.22 

1.20 

3.3x1 0 6 



0 

1.03 


1.05 

1.05 





15 

1.04 


.993 

.999 





75 

1.52 


1.23 

1.16 


! 



90 

1.56 


1.33 

1.33 

6.6xl0 6 




0 

1.05 


1,05 

1.05 





15 

1.11 


1.07 

1.04 





75 

1.73 


1.41 

1.36 

> 

' 



90 

1.77 


1.57 

1.57 

8.2 

O 

X 




0 

1.11 


1.11 

1.11 





15 

1.20 


1.13 

1.13 





75 

1.85 


1.57 

1.53 


' 

1 


90 

1.90 


1.75 

1.69 









PANEL -6 

Re / M 

00 

e ( cm) 

M 

00 

<f> 

(DEG) 

P /P 
max 7 fp 

3.3xl0 6 


79 

2. 

.4 


0 


1.02 








15 


1.07 








30 


1.14 








60 


1.27 

1 







90 


1.33 

10,xl0 6 

1 

I 





0 


1.01 








15 


1.08 








30 


1.18 





i 



60 


1.32 





1 



90 


1.40 





4.5 


0 


.993 








15 


1.05 








30 


1.16 



\ 


1 



60 


1.47 

1 







90 


1.58 
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TABLE 6-V AVERAGE WAVE PRESSURE 


PANEL -1 PANEL -2 

(e = 1.06cm) (e = .525cm) 


Re 

CO 

M 

00 


$ (DEG) 

V p f P 

P av //p fp 

1.3xl0 6 

10 

.3 

0 

.874 

.965 





15 

.929 

.928 





75 

1.00 

.951 


' 



90 

.959 

.972 

3. 3x1 0 6 



0 

.939 

.931 





15 

.904 

.922 





75 

1.07 

.978 

' 

' 



90 

1.05 

1.03 

6.6xl0 6 



0 

.982 

.969 





15 

.953 

.971 





75 

1.18 

1.06 

\ 

( 



90 

1.16 

1.14 

8. 2xl0 6 



0 

1.04 

1.04 





15 

1.02 

1.03 



\ 


75 

1.27 

1.15 

'1 

f 

1 


90 

1.25 

1.25 


RejM e(cm) 


<f> (DEG) 


PANEL -( 
P av /P fp 


3.3xl0 6 .79 2.4 


0 

15 

30 

60 

90 


1.01 

1.03 

1.04 
1.06 
1.07 


lO.xlO 6 


0 

15 

30 

60 

90 


1.01 

1.02 

1.04 

1.07 

1.08 


4.5 


Y 


0 

15 

30 

60 

90 


.987 

1.00 

1.02 

1.11 

1.13 


PANEL -3 
e = .405cm) 

av^fp 


.926 

.962 

1.04 

1.01 

.944 

.911 

.993 

1.09 

.984 

.976 

1.08 

1.19 

1.05 

1.06 
1.18 
1.29 
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panel -1. The results of the calculations are presented in Figures 6-1 and 6-2. 

For both cases it was found that the maximum surface pressure could be predicted. 

For the supersonic case, a value of P m=v /Pr D = 1.36 is predicted compared to 
measured ratios of 1.34 and 1.4 for Re^/M = 3.3 X 10 b and 10 X 10 , respectively . 

For hypersonic flow a maximum value of 1.57 is predicted compared to an experimental 
value of 1.60 for Re /M =3.3 X 10 6 . However, the flow fields required for the pre- 

oo 

diction of the pressures are inconsistent with the measured results. For M^ = 2.4, 
the flow at the top of the wave is required to expand through an angle of 35° over 
the wave to a static pressure value into the shock of P-j/Ppp ~ -125, and for = 10.3 
the required static pressure value is .101 for an expansion angle of 33°. The 
measured separation angles are 5.2° for = 2.4 and 28° for 10.3. Hence, the 
necessary separation angles are much less than detected by experiment. In addition 
the extremely low static pressures required for the flow into the shock wave are 
much less than the measured values which were on the order of P-|/Ppp = 1. Hence, 
the peak pressures measured on the compression surface of the waves cannot be 
explained on the basis of the local internal flow passing through a shock wave. 

An interesting feature of the flow field over shallow waves in thick turbulent 
boundary layers is that the incoming flow adjusts almost instantaneously from wave 
to wave. This was first detected in Ref. 8 in which the heating distributions 
were found to be essentially independent of wave cycle. In order to further in- 
vestigate this phenomenon, a total pressure rake was placed behing the last wave 
of panel -1, for 0 = 90° as shown in Figure 6-3. Total pressures obtained from 
the rake data are compared with maximum pressures measured on a wave in the center 
of the panel as a function of Reynolds number in this figure. The total pressures 
from the rake data were obtained at the measured reattachment height on the wave. 

As shown in Figure 6-3, the two pressures agree quite well which indicates that the 
mixing process over each wave is almost identical, and that the maximum pressure 
on the wave is approximately equal to the total pressure of the streamline that 
stagnates on the compression surface. 

Based on these results, it was decided to determine if the maximum surface 

wave pressure could be related to the total pressure associated with the streamline 

at the wave height. It was found that the ratio of the maximum surface pressure 

to the total pressure of the streamline at the wave height is a function of the 

Mach number of the streamline at the wave height, as shown in Figure 6-4. Hence, 

the ratio P /p t can be considered a total pressure loss associated with the 
max 
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EXPERIMENT 

=1.40 (R e /M=10X10 6 ) 


FIGURE 6-1 IMBEDDED SHOCK FLOW MODEL FOR 

PANEL -6, M = 2.4, 0 = 90° 

00 



NORMAL SHOCK 



FIGURE 6-2 IMBEDDED SHOCK FLOW MODEL FOR 

PANEL -1, M = 10.3, R e /M = 3.3X10 6 , 0 = 90° 

00 t *oo 



*» j * x 10-6 (e t/pfp)* w p f P 
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FIGURE 6-3 


COMPARISON OF MEASURED MAXIMUM PRESSURE AND PITOT PRESSURE 
BEHIND PANEL AT REATTACHMENT HEIGHT, PANEL -1,0= 90° 







turbulent mixing process over the wave. It is not surprising that the pressure 
ratio is such a strong function of Mach number. The Mach number is known to 
strongly affect the total pressure loss in turbulent jets and wakes. 

Shift in Peak Heating and Pressure 

The panels considered in the present study were instrumented in much more 
detail than those previously tested. As a result the measured heating and pressure 
distributions on panels -1, -2, -3, and -6 are much more exact than those on panels 
-4 and -5. Due to the more detailed instrumentation, certain features of the data 
have been detected that were not noticeable in the previous study. For example, 
consider the peak heating and peak pressure locations for 0 = 90° on panels -1, -2, 

-3, and -6 as shown in Table 6- V I . For panels -2 and -3 there is a large shift in 
the peak pressure location. These data are for hypersonic flow at the lowest 
Reynolds number tested. Table 6- V I shows that several flow parameters have 
distinct changes for these conditions which could result in this shift. The shift 
occurs when the Mach number at the wave height becomes subsonic. For supersonic 
flow at the wave height a separation shock could occur just downstream of the wave 
peak. For subsonic flow at the wave height the flow could expand further in the 
valley before separation occurred. In addition to the Mach number change, the sub- 
layer thickness becomes greater than the wave height, and the hydraulic diameter, 
dft*, of the wave shape becomes less than two times the laminar sublayer thickness 
when there is a shift in the peak pressure location. The changes in the later two 
ratios could affect the manner in which the flow field develops over a wave. The 
results presented in Table 6-VI show that a definite change in the flow field over 
the waves has been detected as a function of changes in the local flow conditions. 

The mechanism controlling these changes needs further investigation. It should 
not be ruled out that nonuniformity in the wall boundary layer flow as a function 
of Reynolds number could account for some of these changes. The shift occurred 
only for pressure, and the peak heating location was not affected by the afore- 
mentioned changes in the flow conditions. 

Large shifts in peak pressure locations were also detected for other flow 
angles for the low Reynolds number data, as shown in Figures 6-5 and 6-6 for panel 
-1 for 0 = 15° and 75° as a function of Reynolds number. Hence, there is a 
definite shift in peak pressure location even for supersonic flow at the wave height 
for certain wave shapes. Similar results have been detected by Howell (Ref. 14) in 
his study of flow over notches. Howell found that the separation point for notch 
flow was both a function of the notch shape and the local edge Mach number. 

*The hydraulic diameter, d. , is described in Table 6- V 1 1 . ORIGINAL PAGE IS 
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PANEL 

M 

CO 


Re„/ M Xl(f 6 

6 S (cm) 
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£ 
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e 

d h 

•26? 
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10. 

.3 

1.3 
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1.44 
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1.31 

.8 

17 

.8 
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3.3 

.381 

1.93 

.36 

2.53 
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.279 
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.26 

3.45 
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.19 

4.74 
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1 

> 

2 



1.3 
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.81 

1.41 
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.75 
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.381 

1.20 
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' 
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.203 
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.82 

> 

f 
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1.3 

.737 

.73 

1.82 

.52 

.73 
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3.3 

.381 

1 .00 

.94 

1.01 

.92 

.93 





6.6 

.279 

1.27 

.69 

1.38 

.91 

.91 


f 

) 


8.2 

.203 

1.43 

.50 

1.89 

.90 

.90 

6 

2. 

4 

3.3 

.028 

1.5 

.035 

7.89 

.75 

.81 

1 

2. 

,4 

10.0 

.013 

1.5 

.017 

17.0 

.75 

.80 


4. 

,5 

10.0 

.062 

1 .95 

.079 

3.57 

.76 
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TABLE 6- V I PARAMETERS AFFECTING PEAK HEATING AND PRESSURE LOCATIONS , 0 = 90° 
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TABLE 6-VII DEFINITION OF HYDRAULIC DIAMETER 
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FIGURE 6-5 SHIFT IN PEAK PRESSURE LOCATION, M =10.3 

00 
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FIGURE 6-6 SHIFT IN PEAK PRESS 


( 





Howell's correlation of the separation point is shown in Figure 6-7 in which the 
data are given as a function of the Mach number and notch angle. For a given notch 
angle and subsonic flow, the separation point moves upstream as the Mach number 
increases. - The separation distance reaches a minimum for sonic flow, and then 
increases as the Mach number becomes supersonic. For a given local edge Mach 
number, the separation distance decreases as the notch angle is increased. 

A similar plot of the present data for corrugated shapes is shown in Figure 
6-8. In plotting these data, the separation point on the expansion surface was 
assumed to be the mirror image of the measured reattachment point (point of surface 
peak pressure) on the compression surface. Judgement was used in constructing the 
faired distributions due to the lack of data in the subsonic and sonic flow range. 
Nonetheless, these distributions seem reasonable considering Howell's measured 
distributions. Hence, it appears that the Mach number at the wave height controls 
the separation point location on a corrugated shape in a thick turbulent boundary 
layer. For a given Mach number at the wave height, the separation distance de- 


creases as the parameter 


rL 

c-2 


is decreased. 


Correlation of Peak" Heating 

The previous, study. Refs. (8) and (9), served as a guide in choosing the 
parameters which would be useful in developing the peak heating correlations. An 
automated multiple regression technique designed to fit multiple variables was 
employed to help determine the importance of the input parameters and obtain con- 
sistent nonbiased correlating equations. The computer program, which is described 
in detail in Ref. 15, is a stepwise Multiple Regression Analysis (MRA) which pro- 
vides information as to the adequacy of candidate correlation functions and the 
equation coefficients. The MRA computes a series of multiple linear regression 
equations in a stepwise manner. At each step, one parameter is added to the 
equation. The variable that is added is the one which makes the greatest reduction 
in the variance about the mean. 

Formulation of the first peak heating correlation consisted in developing a 
good analytical fit to all the thick turbulent boundary layer data obtained on 
panels -1, -2, -3, -4, -5, and -6. These maximum heating data ranged from a value 
of 1.37 to 4.63. This range includes all the peak heating data except for corrugati 
angles of 0°. Several good fits were formulated. The best fit of the data has a 
standard error of estimate of .324 and is given by the equation: 
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FIGURE 6-7 SELF ADJUSTMENT OF THE ' SEPARATION POINT FOR THIN TURBULENT 
BOUNDARY LAYER FLOW OVER NOTCH 


rj^tsioiao 



• J 1*0 1.5 m 2.0 

e 

FIGURE 6-8 SELF ADJUSTMENT OF THE SEPARATION POINT FOR THICK TURBULENT BOUNDARY LAYER 
FLOW OVER WAVY WALL 

A 


2.5 


0 ) 


h max . r , . v 390 M 1 . 077 * 329 

-1 + C 1 (sincx r ) M £ R 6££ 




.196 


< sin W 


.325 


where = e~' 677 

During the course of developing this correlation equation, it was discovered 
that several of the terms could be eliminated without greatly reducing the accuracy 
of the resulting correlation. Removing the displacement thickness, momentum thick- 
ness, and the maximum wave surface angle, the following correlation equation was 
obtained for the peak heating for all the thick turbulent boundary layer data: 


h 

max 

h FP 


1 + C , 


(sina r ) 


.240 M 1.316 



.231 


( 2 ) 


951 

where C2 = e’ . This simpler equation has a standard error of estimate of .350. 

Equation (1) is compared with all the thick turbulent boundary layer data in 
Figure 6-9. Equation (2) is not graphically compared with the data. However, the 
agreement would be comparable since the standard error of estimate is about the 
same for both equations. 

Prior analysis (Ref. 8) using the MRA solution for peak heating showed that 
good data fits for thick turbulent boundary layer data were in poor agreement with 
data sets for thin turbulent boundary layers. As an extension of the previous 
correlation activity, all the thin boundary layer data were incorporated into the 
data bank containing the thick boundary layer data. The MRA analysis was then 
applied to the resulting data bank containing all the thick and thin boundary layer 
peak heating data. A successful correlation equation was then derived for a range 
in data from 1.37 to 7.6. This equation is: 



1 + C 3 (sina r ) M £ 1 - 107 





.344 


(3) 
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FIGURE 6-9 CORRELATION OF MAXIMUM HEATING FOR ALL 
THICK BOUNDARY LAYER DATA 




, r 2.813 

where C^ = e 

This correlation, which has a standard error of estimate equal to .557, is compared 
with all the data in Figure 6-10. The agreement between the measured and predicted 
values is good except for a few thin boundary layer data points at the highest 
heating conditions. 

As with the analysis for all the thick boundary layer data, it was found that 
several correlating parameters could be eliminated from equation (3) in order to 
simplify the correlation without greatly reducing its accuracy. A simpler 
correlation form for all the peak heating data is: 


h 

max 

h FP 


1 + C. (sina ) ,348 M 1,397 
4 r e 



(!) ' 678 


(4) 


where = e 4 ' 131 ^ 

The standard error of estimate for this correlation is .579. Agreement between 
the correlation predictions and the data is comparable to that shown in Figure 6-10 
because the standard errors of estimate for equations (3) and (4) are very similar. 
Correlation of Average Heating 

An important factor which must be considered in the study of corrugated panels 

is the effect of surface roughness on average heating. Correlations and analyses 

for the prediction of average heating to rough surfaces have previously been 

developed for nosetip performance during reentry conditions. A correlation 

developed by Powars (Ref. 16) has been in widespread use for the prediction of 

surface roughness heating. Data obtained on all six panels of the present study 

are compared with Powars correlation in Figure 6-11. The equivalent sand grain 

roughness, k g , used in this correlation was calculated using Dirling's analysis 

(Ref. 17). The data do not follow the trend of the correlation for any range of 

the independent variable. In general the data fall well below the correlation line 

for 10k / 6 >10 and have considerable scatter about the correlation line for 

e s 

10k / 6 <10. 

e s 

Another correlation for average heating to rough surfaces was developed by 
Grabow and White (Ref. 18) using the same heating data upon which the Powars 
correlation is based. Data from the present study are compared with the Garbow 
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FIGURE 6-11 AVERAGE HEATING COMPARED WITH POWARS CORRELATION 
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and White correlation in Figure 6-12. In this correlation the average heating is 
plotted versus k e /6*. This correlation works quite well for k e /6* greater than 
about 10. The present data fall in the range for ke/6* less than about 1. These 
data represent the range for which the roughness height is much less than the dis- 
placement thickness. The correlating parameter k e /6* does not represent an 
adequate variable which can be used in predicting the average heating for very 
thick turbulent boundary layers. 

Based on the comparisons shown in Figures 6-11 and 6-12, it was concluded 
that no reliable correlation existed which could be used to represent average heat- 
ing to rough surfaces in thick turbulent boundary layers. Some success had pre- 
viously been obtained in Ref. 8 using Dirling's technique to predict the average 
heating for the flow normal to the corrugations. However, there was an uncertainty 
as to which geometric height (half-wave height or total wave height) must be used 
in order for Dirling's method to match the experimental results. The reason for the 
uncertainty was not understood and needed more investigation. It was felt that a 
detailed analysis of the present data using Dirling's approach was beyond the scope 
of this study. Rather, it was decided to focus effort toward developing a correlati 
equation using the MRA approach. 

It was found that an excellent correlation of the average heating data could 
be obtained using the corresponding maximum heating data along with several flow 
and geometric parameters. This correlation is: 



, .454 .174 0 .012 x * .005 

'* (f) (r) (f) 


where C 5 = e" 2 ‘^ 74 . The measured data for this correlation range from .836 to 1.85. 
The correlation has a standard error of estimate of .081. Equation (5) is compared 
with the present data in Figure 6-13 for all corrugation angles except 0 = 0°. The 
agreement between the correlation and the data is quite good. 

Further analysis showed that the complexity of equation (5) could be greatly 
reduced by eliminating several correlating parameters without severely penalizing 
the accuracy of the curve-fit to the data. The best simplified correlating form 
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FIGURE 6-12 AVERAGE HEATING COMPARED WITH GRABOW AND WHITE CORRELATION 
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FIGURE 6-13 CORRELATION OF AVERAGE HEATING FOR ALL 
THICK BOUNDARY LAYER DATA 
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for the average heating to rough surfaces in very thick turbulent boundary layers 
is : 



__ -I oO£ 

where Cg = e . This simplified correlation has a standard error of estimate 
of .088. The agreement between this correlation and the data is essentially as 
good as that shown in Figure 6-13. 

Geometry of Separation Region 

Inherent in the maximum heating correlations is the assumption that the geometry 
of the separated flow field in the valley of a wave is known so that the local 
wave surface angle can be computed at the point of maximum heating. Prior to the 
present study, the local wave surface angle at maximum heating was assumed to be 
known from experiment. If the angle was not known, it was suggested that the 
maximum surface angle on the wave be used in order to obtain a conservative 
solution for the maximum heating. 

In order to improve the usefulness of the developed correlations, the experiment- 
al data for the surface angle at the maximum heating point have been correlated for 
all the thick boundary layer data obtained on panels -1, -2, -3, -4, -5, and -6 in 
supersonic and hypersonic flow. The data for thin boundary layers were not in- 
cluded in the geometry correlations since the local angles could not be accurately 
computed for the data that had been published in the literature. The following 
equation was found to give a good fit to the thick boundary layer geometry data: 


s i na r = 



.320 


(7) 


6 026 

where Cy = e°* . This correlation has a standard error of estimate equal to .073. 

Again, further examination showed that the equation could be greatly simplified 
without reducing the accuracy of the results. Upon simplification, the following 
equation resulted which represents a good fit of the data: 
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/ V 1.373 .097 

Sina r C 8 (U J Ree e (8) 

/ . x .63 
(s1n “max> 

where Co = e *^ . This correlation has a standard error of estimate of .077. 
Equation (7) is compared with the data in Figure 6-14. Except for a few points, 
the agreement between the measured and predicted results is fairly good. About 
the same agreement would exist between the data and equation (8). 

Variation of Local Heating with Pressure 

In many engineering studies, a heating-pressure relationship is used to pre- 
dict the maximum heating. The variation in the heating rate with the local pressure 
measured in this study is shown in Figure 6-15. All the thick boundary layer data 
are shown in this figure for 0 = 90°. Two separate relationships exist between 
pressure and heating. In the separated region, the heating increases while the 
pressure is fairly constant. In the vicinity of the attached flow region, the 
heating ratio is approximately proportional to the square of the pressure ratio. 

The reattachment heating for thin turbulent boundary layers varys as the pressure 
ratio to about the eight- tenth power. Hence, in comparison to thin boundary 
layers, the present results for thick turbulent boundary layers over corrugated 
surfaces show that the increase in heating over the corresponding flat plate value 
is much greater than the increase in the pressure over the flat plate value. 
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7.0 THERMAL FLIGHT ASSESSMENT 

Using the data and correlations described in the previous sections, this 
section examines the feasibility of applying surface corrugation stiffened panels 
to advanced space transportation system (STS) vehicles without serious risk to flight 
performance. Factors considered in the evaluation are the impact of corrugations 
on thermal response and their impact on aerodynamic performance, both within con- 
straints established by structural dynamic criteria and materials technology. 
Optimization of corrugation shapes with respect to cost, weight or system per- 
formance was not considered. 

Objectives of the thermal assessment were to determine, from application of 
the MRA correlations for peak and average corrugation heating, the peak TPS skin 
temperatures to be expected during reentry and the insulation thickness required 
to maintain a 177°C (350°F) limit for the primary structure temperature excursion 
at several representative points on the SSTO vehicle. From these results it was 
the further objective to gauge, by the predicted temperature and insulation thick- 
ness departures from those obtained for the corresponding smooth body case, where 
on the SSTO vehicle the corrugated panel TPS can be successfully applied without 
excessive penalties in insulation thickness or skin material requirements. 

The peak corrugation temperatures were calculated by applying to the corrugated 
skin an assumed spatial heating distribution scaled to result in the same peak 
heating value as that predicted by the MRA correlations. This calculation was 
repeated at each location on the vehicle for two different materials (depending 
upon the expected temperature range) and three skin thicknesses to determine the 
effect of variations in lateral skin conductance upon the predicted peak temperature. 

Insulation thicknesses were determined at each location on the vehicle for 
both the smooth body and the corrugated skin cases using, in the latter case, the 

correlated value for average corrugation heating to augment the smooth body value. 

3 

This calculation was performed for only one skin thickness and assumed 56.1 kg/m 

O 

(3.5 Ibm/ft ) Microquartz insulation. 

Throughout this analysis, the assumed corrugation shape is that shown in Figs. 
7-5 and 7-6 which is the configuration recommended by the strength and structural 
dynamics considerations outlined in Section 9.0. 

SSTO Vehicle and Trajectory 

While the study results are intended to be applicable to advanced STS vehicles 
in general, to facilitate the analysis a specific SSTO (Single-Stage-To-Orbit) 
vehicle was chosen as the focal point for the flight evaluation. The characteristics 
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of this vehicle were developed in Ref. (19) and shown in Figure 7-1. The SSTO is 

designed to be sled-launched and to land horizontally. Its length, 63m (206 ft), 

and wing span, 40m (131 ft), are roughly twice the corresponding dimensions of 

2 

Shuttle. The reentry planform loading of SSTO is approximately 129 kg/m 
(26 lb/ft^) compared with 264 kg/m^ (55 lb/ft^) for Shuttle. 

The trajectory flown by the SSTO vehicle from Ref. (19) is shown in Figure 
7-2. The angle of attack is 50° at the time of deorbit, decreasing to 30° 2300 
seconds later and remaining at that value for 2700 seconds. The period of sign- 
ificant aerodynamic heating and peak loads occurs within the period during which 
the SSTO angle of attack is 30°. The fact that the angle of attack remains 
constant during this period results in simplification of the corrugated panel 
design evaluation. 

Vehicle Locations Selected 

Figure 7-3 illustrates the location of the points of interest on the SSTO for 
which the flow field and its effects upon corrugation heating were investigated. 
These locations were selected to be representative of areas in which it is expected 
that the corrugated panel TPS might be applicable. In addition, the points were 
chosen to reflect a wide range of smooth wall heating and boundary layer parameters 
in order to define rough limits for the application of corrugated panels on the 
vehicle and for the correlations developed to predict peak and average corrugation 
heating. An important consequence of the choice of points for the flow field 
evaluation is the wide variation in local flow angle relative to the vehicle axis, 
and thus the possible variation of the flow direction relative to the assumed axis 
of the TPS corrugations. As was previously noted in the data correlation, this 
corrugation "angle of attack" is one of the more important parameters affecting 
peak heat transfer to the corrugation. 

Wind tunnel oil flow results for Shuttle orbiter type configurations at 30° 
angle of attack (Figure 7-4) were used as a basis for determining the local flow 
angle relative to the vehicle roll axis. These suggest that flow angles for the 
SSTO geometry are probably less than 10° near the lower body centerline, less than 
20° at the lower body off-center point, and 20° on the lower wing. While the 
flow on the upper surfaces is predominantly separated and shows evidence of vortex 
impingment, the local flow angles appear to be no greater than 60° at any point 
on the upper body, upper wing or side body forward of the wing. 
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FIGURE 7-1 SSTO VEHICLE CONFIGURATION ASSUMED FOR THE PRESENT ANALYSIS 
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FIGURE 7-3 VEHICLE LOCATIONS CONSIDERED FOR APPLICATION OF THE CORRUGATED PANEL TPS 
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FIGURE 7-4 OIL FLOWS USED TO ESTIMATE LOCAL FLOW ANGLES RELATIVE TO THE VEHICLE ROLL AXIS 
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FIGURE .7-4 CONTINUED LEEWARD SURFACES 









The MRA correlations for peak heating and associated estimates of the spatial 
distribution of seating and temperature on the corrugation were evaluated for these 
assumed maximum corrugation angles to determine the peak temperature penalty in 
regions for which this angle is poorly defined, i.e., for the high local flow 
angles typical of leeward surfaces. Alternately, the possibility of a significant 
reduction in peak corrugation temperature was allowed for regions in which close 
alignment of the corrugation with the local flow may be realistic (typical of wind- 
ward surfaces). 

Evaluation of the MRA average heating correlation and estimates of insulation 
thickness were also confined to the estimated maximum local flow angles that would 
result from simple alignment of the corrugations with the vehicle roll axis, i.e., 
60° on leeward and side surfaces, 10° on the lower body centerline and 20° for 
the lower body off center point and lower wing. 

Inviscid Flow Models 

Trajectory time dependent boundary layer edge conditions and smooth body heat 
transfer coefficients for each point of interest on the vehicle were calculated 
using the Miniver version of the JA-70 aerodynamic heating code (Ref. 20). 

On windward surfaces, it was assumed that the local pressure is approximated 
by that predicted for cone flow. In regions of separated flow on the leeward 
side, turbulent attached flow was assumed for the purposes of modeling and a 
turbulent design factor (heat transfer coefficient multiplier) applied to obtain 
smooth body temperatures that are in agreement with those determined experimentally 
(Ref. 21). 

The flow models constructed and additional assumptions made are as follows: 
Lower body 

o shock wave angle is obtained from a table of sharp cone shock angles 
as a function of upstream Mach number and a cone half angle equal to the 
time dependent angle of attack of the vehicle. 

o local pressure is determined from a table of sharp cone real gas pressure 
coefficients as a function of upstream Mach number and cone half angle. 

o local smooth body enthalpy-based heat transfer coefficient is calculated 
using a real gas extension of the Spalding Chi flat plate method. 

o correction to account for streamline divergence is made to the geometric 
running length by using the real gas crossflow technique developed by 
Baranowski (Ref. 22). 
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Lower wing 

o flow model is the same as that used for the lower body except that the 
geometric running lengths are shorter and the crossflow correction is 
omitted. 

Upper body and side body 

o shock wave angle is obtained in the same manner as that used for the lower 
surfaces with the cone half-angle again assumed to be equal to the time 
dependent angle of attack of the vehicle. 

o local pressure is determined by specifying the local pressure coefficient 
C p - 0. 

o local heat transfer is calculated using the Spal ding-Chi flat plate 
method with the turbulent design factor set equal to 2.0. 
upper wing 

o shock wave angle is assumed equal to the effective sweep angle of the 
wing leading edge. 

o local pressure coefficient Cp = 0. 

o local heat transfer is calculated using the Spalding and Chi method with 
the turbulent design factor set equal to 3.5 for the upper wing forward 
point and 2.0 for the aft point. 

Boundary Layer Model 

Local edge conditions calculated by Miniver were input to a simple 1/7 power 
law flat plate boundary layer model based on the work of Walker and Schumann 
(Ref. 23) to obtain the bulk and internal boundary layer parameters required to 
evaluate the correlations for peak and average corrugation heating. The 
following relations were used: 

(ratio of boundary layer thickness to momentum thickness) 

f = (n+1) + + l) (1+ArV) 


n 

h 


w 


aw 


A 

r 

V 


= velocity power law exponent 
= enthalpy at the wall 
- h * r u » 2 

CO ” 

2g c J 

= 1.69 (a constant) 

= recovery factor 


2g c J h M 
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(ratio of displacement thickness to momentum thickness) 

f ■ - 1 + ( ii r L r~ + ^ ) o+Brv) 

' aw ' 

B = 1.16 (a constant) 

(momentum thickness) 



k = 0.037 for flat plates 

The resulting velocity profile and the Crocco expression for the enthalpy- 
velocity relationship through the boundary layer 

h - h „ + < h a W -\,> < 5 - ) - < h a W - h J <£- > 2 

00 CO 

were used to determine internal boundary layer parameters (R £ , M £ , etc.) at the 
corrugation height. 

The laminar sublayer thickness was evaluated from Reda's expression (Ref. 24). 

6 = 11y w a w 

s Y w P w M oo'4 C f/2 

where the skin friction coefficient is calculated in Miniver using Reynolds 
analogy. 

A summary of the flow conditions calculated using these inviscid and boundary 
layer models for each of the ten vehicle locations investigated is shown in 
Table 7-1. 

Numerical Boundary Layer Solution 

While the approximate solution used in this study to obtain bulk and internal 
boundary layer parameters is easily evaluated and thus convenient for repetitive 
calculations, it represents a great simplification of the actual viscous flow. 

For this reason, comparison was made with a numerical solution developed by Fivel 
(Ref. 25). This is an "exact" equilibrium finite difference solution to the com- 
pressible turbulent boundary layer equations for a real gas and has previously 
been applied (Ref. 26) to a 30° cone for freestream conditions that approximate 
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TABLE 7-1 SUMMARY OF PREDICTED FLOW CONDITIONS ON SSTO 



those of the SSTO reentry trajectory near peak heating (Table 7- 1 1 ) . Results of 
the cone boundary layer solution, at a point chosen to correspond to the SSTO 
lower body forward location, were used to evaluate the MRA correlations for peak 
and average corrugation heating, and comparison made with the values obtained using 
the approximate solution. As seen in Table 7- II agreement is quite good in spite 
of the sizable difference in displacement thickness predicted by the two methods. 
Scatter about the mean of the four peak values is only +13%, -3% and the greatest 
departure of the peak value based on the approximate boundary layer from that 
based on the exact solution is -14%. 

The discrepancy in displacement thickness is, at least in part, attributable 
to the correction made to the turbulent running length via the Baranowski cross- 
flow technique (Ref. 22). This adjustment was included in the approximate cal- 
culations to account for the spanwise flow on the SSTO forward lower surface and 
the consequent deviation of the boundary layer development from that expected for 
a flat plate. It would be expected that an "exact" solution based on a geometry 
that corresponds to the SSTO configuration would result in an increase in the pre- 
dicted displacement thickness over the conical prediction and a consequent decrease 
in the peak corrugation heating to a value nearer that obtained from the approximate 
calculations. Of course, the mismatch in free stream conditions may also con- 
tribute to the observed differences. 

The importance of this comparison is two fold. First, it demonstrates that 
the approximate analysis used throughout the study results in reasonable boundary 
layer parameters at least on the windward surfaces and, second, that evaluation of 
the MRA correlations using these parameters yields values of peak corrugation 
heating that are in agreement with those predicted using the more rigorous 
numerical analysis. 

TPS Skin Lateral Conduction Model 

Since the correlations developed during the present study predict only the 


peak heating on the windward^ side of the corrugation and do not address the spatial 
distribution of heating ovijf its wavelength, it is necessary, if advantage is to 
be taken of lateral conduction in the skin, to make some assumption about the "shape" 
of this distribution. This done, the actual peak temperature of the surface may 
be calculated. Further, it was desirable to perform these calculations for a 
range of skin thickness to determine the importance of skin lateral conduction 
in reducing the peak temperature. 
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FREE STREAM CONDITIONS 
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A review of the experimental heating profiles disclosed that they appear to 
be of two genre, one having a broad region of relatively high heating rate near 
the peak and only narrow regions of heating at or below the flat plate value. 

This distribution is typical of the -6 type corrugation (Ref. 10) shown in Figure 
3-1. Corrugations of a more nearly "sinusoidal" cross-section such as those of 
the -5 type show a relatively narrow region of high heating rates near the peak 
and broad areas in which the heating is significantly below the smooth body value. 

Since it is not entirely clear that this implied correspondence between the 
two types of corrugations and their respective heating profiles is real, profiles 
representative of both types were taken as reference values from which to scale 
heating distributions corresponding to the predicted peak values. These were then 
used to calculate two values for the peak temperature expected on the corrugation 
(lateral conduction effects included), and thus estimate the magnitude of the 
effect of uncertainties in the heating distribution. 

The peak temperatures were calculated using a 20 node, one dimensional thermal 
model which includes conduction laterally through the metallic TPS skin, aerodynamic 
heating at each node based upon an evaluation of the MRA peak heating correlation 
at flow conditions predicted by Miniver for trajectory peak heating, and reradiation 
from the skin outer surface (emissivity = 0.8) to a 3 K radiation sink. Figures 
7-5 and 7-6 show a schematic of the model, typical applied heating distributions, 
and the resulting temperature profiles for three different skin thicknesses. 

These calculations assume equilibrium conditions and an adiabatic back surface; 
however, during the course of the transient heating calculations using Miniver, it 
was noted that the TPS skin does, in fact, very nearly attain the radiation 
equilibrium temperature even for the relatively large skin thicknesses assumed. 

For Miniver calculations on the lower body and lower wing a skin thickness of 
0.254 cm (0.1 in) was assumed for the purposes of computing the inviscid flow 
field since this caused a departure of the wall temperature from the radiation 
equilibrium temperature of only a few degrees and allowed calculations to be made 
at less frequent (more economical) time intervals. The lower heating rates 
characteristic of the vehicle upper surface required a thinner assumed skin 
(0.076 cm) to maintain the skin temperature at a value near the radiation 
equilibrium temperature and a corresponding reduction in the calculation interval 
to insure a stable solution. In either case, for the present purposes of skin 
material selection, it may be assumed that the typically much thinner (0.041 cm) 
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LOCAL CORRUGATION TEMP. (K) 
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PROFILE "A" - HEATING DISTRIBUTION 



FIGURE 7-5 PREDICTED CORRUGATION TEMPERATURES FOR LOWER 

FUSELAGE FORWARD LOCATION (HEATING PROFILE "A") 
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LOCAL CORRUGATION TEMP. (K) 







TPS skin to be used on a flight vehicle is in thermal equilibrium at the radiation 
equilibrium temperature. As a result, it may further be assumed that the 
equilibrium lateral conduction model accurately reflects (within the uncertainty 
imposed by the assumed heating profile) the actual peak temperatures to be expected 
during transient heating of the TPS. 

Results of these calculations at each point of interest on the SSTO for both 
heating profiles, skin thicknesses of 0.0, 0.015, and 0.041 cm (0.0, 0.006, and 
0.016 in), and several different materials depending upon the temperature range 
considered are shown in Table 7- 1 1 1 . The correlated ratios of peak and average 
heating and the corresponding smooth body heat transfer coefficients are listed in 
Table 7-IV. 

The effects of radiation and conduction heat transfer from the backside of 
the TPS skin to the insulating medium will be discussed in the following section. 
TPS Insulation Sizing 

To assess the impact of the increased average convective heating to the 
corrugation surface upon. TPS insulation requirements, a one-dimensional transient 
conduction model was constructed based upon the following assumptions: 

o a fixed TPS outer skin thickness of 0.041 cm (0.016 in). 

3 3 

o 56.1kg/m (3.5 Ibm/ft ) Microquartz insulation (20 nodes) with temperature 

and pressure dependent conductivity. 

o 0.254 cm (0.1 in) 2024-T4 aluminum backside to simulate the vehicle 
primary structure. 

o a primary structure transient temperature limit of 450K (350°F). 

o an initial temperature of 255K (0°F). 

o adiabatic conditions at the primary structure back surface. 

Smooth body transient heating rates were read into the model directly from a 
Miniver output tape consisting of time dependent recovery temperature-heat 
transfer coefficient pairs that result in the same heating rate as the correspond- 
ing recovery enthalpy and real gas heat transfer coefficient. The insulation 
thickness was then varied in trial and error fashion to determine the smooth body 
insulation thickness required to limit the primary structure temperature 
excursion to 450K. Corrugation average heating rate increments determined from 
the MRA correlations were then input as multipliers on the smooth body values and 
the trial and error insulation sizing repeated to determine a new thickness 
required to maintain the 450K limit. 
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I 

rv> 

o 



A 


CORRUGATION TEMPERATURE (K) 




PROFILE 

A 

PROFILE B 


Z(cm) 

MAX 

MIN 

MAX 

MIN 

LOWER FUSELAGE FWD 0 = 10° 

.041 

1405. 

1278. 

1377. 

1222. 


.015 

1431. 

1255. 

1414. 

1196. 

(SMOOTH BODY TEQ = 1257 K) 

.000 

1466. 

1212. 

1466. 

1169. 

AFT 0 = 10° 

.041 

1312. 

1237. 

1293. 

1206. 


.015 

1328. 

1224. 

1317. 

1191. 

(SMOOTH BODY TEQ = 1224K) 

.000 

1355. 

1199. 

1353. 

1177. 

OFF CENTER 0 = 20° 

.041 

1337. 

1242. 

1314. 

1199. 


.015 

1358. 

1224. 

1344. 

1169. 

(SMOOTH BODY TEQ = 1224K) 

.000 

1389. 

1191. 

1388. 

1161. 

LOWER WING FWD 0 = 20° 

.041 






.015 


NOT 



(SMOOTH BODY TEQ = 1324K) 

.000 


EVALUATED 



AFT 0 = 20° 

.041 

1506. 

1303. 

1466. 

1207. 


.015 

1541. 

1267. 

1519. 

1159. 

(SMOOTH BODY TEQ = 1271 . K) 

.000 

1584. 

1191. 

1582. 

1107. 


TABLE 7- II I PREDICTED CORRUGATION PEAK AND MINIMUM TEMPERATURE 
FOR WINDWARD SURFACES/TD-Ni20Cr 



LOWER FUSELAGE FWD 0 = 10° 
(SMOOTH BODY TEQ = 1 257. K) 

AFT 0 = 10° 

(SMOOTH BODY TEQ = 1 224. K) 

' OFF CENTER 0 = 20° 


(SMOOTH BODY TEQ = 1224K) 
LOWER WING FWD 0 = 20° 


ft 8 
10 

i! 

F ft 

a ft 


(SMOOTH BODY TEQ = 1324. K) 

AFT 0 = 20° 

(SMOOTH BODY TEQ = 1271 .K) 




TABLE 7 


CORRUGATION TEMPERATURE (K) 



PROFILE A 


PROFILE B 


Z(cm) 

MAX 

MIN 

MAX 

MIN 

.041 

1388. 

1293. 

1354 

1239. 

.015 

1416. 

1268. 

1393. 

1210. 

.000 

1466. 

1212. 

1466. 

1169. 

.041 

1301. 

1247. 

1279. 

1216. 

.015 

1319. 

1232. 

1342. 

1199. 

.000 

1355. 

1199. 

1355. 

1177. 

.041 

1323. 

1253. 

1296. 

1213. 

.015 

1346. 

1234. 

1356. 

1191. 

.000 

1389. 

1191. 

1389. 

1161 . 

.041 

1768. 

1423. 

NOT 


.015 

1828. 

1339. 

EVALUATED 


.000 

1899. 

1070. 



.041 

1483. 

1327. 

1432. 

1238. 

.015 

1522. 

1287. 

1490. 

1187. 

.000 

1584. 

1191. 

1584. 

1107. 


III CONTINUED WINDWARD SURFACES/FS-85 
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CORRUGATION TEMPERATURE (K) 






PROFILE 

_A 

PROFILE 

_B 




Z(cm) 

MAX 

MIN 

MAX 

MIN 

UPPER FUSELAGE 

FWD 

0 = 60° 

.041 

724. 

657. 

682. 

596. 




.015 

749. 

628. 

719. 

556. 

(SMOOTH BODY TEQ = 

609. K) 


.000 

806. 

521. 

806. 

417. 


AFT 


.041 

692. 

638. 

656. 

588. 




.015 

714. 

614. 

686. 

557. 

(SMOOTH BODY TEQ = 

600. K) 


.000 

767. 

528. 

767. 

462. 

UPPER WING 

FWD 


.041 

1080. 

851. 

1024. 

746. 




.015 

1124. 

774. 

1091. 

643. 

(SMOOTH BODY TEQ = 

749. K) 


.000 

1183. 

424. 

1183. 

424. 


AFT 


.041 

795. 

693. 

744. 

615. 




.015 

829. 

651. 

793. 

555. 

(SMOOTH BODY TEQ = 

619. K) 


.000 

892. 

458. 

892. 

348. 

SIDE FUSELAGE 

FWD 


.041 

789. 

693. 

739. 

616. 




.015 

821. 

654. 

786. 

558. 

(SMOOTH BODY TEQ = 

627. K) 


.000 

883. 

489. 

883. 

353. 


TABLE 7 -III CONTINUED. LEEWARD SURFACES/6AL-4V TITANIUM 
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CORRUGATION TEMPERATURE ( K) 






PROFILE 

_A 

PROFILE 

B 




Z(cm) 

MAX 

MIN 

MAX 

MIN 

UPPER FUSELAGE 

FWD 

0 = 60° 

.041 

718. 

663. 

674. 

605. 




.015 

742. 

637. 

708. 

570. 

(SMOOTH BODY TEQ = 

609. K) 


.000 

806. 

521. 

806. 

417. 


AFT 


.041 

687. 

643. 

649. 

594. 




.015 

707. 

621. 

677. 

567. 

(SMOOTH BODY TEQ = 

600. K.) 


.000 

767. 

528. 

767. 

462. 

UPPER WING 

FWD 


.041 - 

1067. 

869. 

1007. 

768. 



o Q 

H" 1 

.015 

1114. 

793. 

1075. 

671 . 

(SMOOTH BODY TEQ = 

749. K) 


.000 

1183. 

424. 

1183. 

424. 



gg 








w £ 







AFT 


.041 

788. 

703. 

733. 

628. 



a ^ 

M O 

.015 

820. 

663. 

779. 

574. 

(SMOOTH BODY TEQ = 

619. K) 

a M 

.000 

892. 

458. 

892. 

348. 









SIDE FUSELAGE 

FWD 


.041 

482. 

702. 

670. 

618. 




.015 

813. 

665. 

696. 

593. 

(SMOOTH BODY TEQ = 

627. K) 


.000 

883. 

489. 

883. 

353. 
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CONCLUDED. 
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LOCATION 

h max /h FP 

h av /h FP 

W v , n 3 

Btu 

o A I U 

cnrK 

Ft 2 -hr-°F 

LOWER FUSELAGE - FWD 

1.97 

1.05 

3.05 

(5.37) 

LOWER FUSELAGE - AFT 

1.56 

.95 

2.72 

(4.79) 

LOWER FUSELAGE - OFF 

1.74 

.96 

2.72 

(4.79) 

LOWER WING - FWD 

5.04 

1.89 

3.83 

(6.75) 

LOWER WING - AFT 

2.65 

1.25 

3.21 

(5.65) 

SIDE FUSELAGE 

5.09 

1.25 

.31 

( .54) 

UPPER FUSELAGE - FWD 

3.84 

1.01 

.27 

( .48) 

UPPER FUSELAGE - AFT 

3.30 

.90 

.25 

( .44) 

UPPER WING - FWD 

9.15 

1.98 

.65 

(1.15) 

UPPER WING - AFT 

5.62 

1.36 

.28 

( .50) 


TABLE 7- IV SUMMARY OF CORRUGATION PEAK AND AVERAGE HEATING RATES AND SMOOTH 
BODY REFERENCE HEATING RATE 
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The correlated average heating multiplier was in all cases assumed to be a 
constant equal throughout the trajectory to the value obtained at the time of 
smooth body peak heating (Table 7- IV) - However, to check the accuracy of this 
approximation, the average heating correlation was evaluated at several points 
in the trajectory prior to and following peak heating for the lower body forward 
location and input to the insulation sizing model as a time dependent multiplier. 

The predicted increase in insulation thickness was approximately 5% which compares 
favorably with 3 % found using the constant multiplier. 

Results of the insulation sizing analysis are shown in Table 7-V. It is seen 
that the greatest percentage increases in insulation thickness occur in areas of 
relatively low smooth body heating and thin boundary layer, e.g., the upper wing 
at 20 % chord. However, on an absolute scale these represent only small increases 
in actual insulation weight since the initial smooth body insulation thickness is 
small. As an example, for the upper wing, use of corrugated panels would result 
in a total insulation penalty on the order of only 70kg (150 Ibm) per wing. 

During the smooth body insulation sizing calculations, it was noted that for 
the vehicle upper surfaces (regions in which the convective heating rate is re- 
latively low) the skin temperature does not attain the radiation equilibrium temper- 
ature. This is obviously a result of heat conduction to the insulating medium, i.e., 
the skin backside is no longer adiabatic. Moreover, it is reasonable to assume 
that a similar effect would occur at the point of peak heating on the corrugation 
if the skin back surface had not been considered adiabatic in the lateral con- 
duction model. A further temperature reduction would be expected to result from 
the transient nature of the heating. The magnitude of these effects upon the 
transient smooth body skin temperature was as great as 34K in one case. In none 
of the cases investigated was the reduction in peak corrugation temperature 
large enough to allow a change to a lower temperature skin material. 

Ideally, of course, the insulation sizing and skin temperature estimate should 
be carried out as part of the same thermal model. This would require knowledge 
of the variation of the spatial distribution of heating on the corrugation during 
the trajectory, a problem not addressed during this study. However, an obvious 
advantage of this approach is that no average heating correlation would be required. 

Three phenomena considered to be of secondary importance have been omitted 
from the insulation sizing analysis: 1) that of "shine through", radiant heat 
transfer caused by the relative transparency of Microquartz (and other silica 
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INSULATION THICKNESS (cm) 
56.1 kg/m 3 MICROQUARTZ 


LOCATION 


SMOOTH BODY 

LOWER FUSELAGE 

FWD 

6.50 


AFT 



OFF CENTER 

5.99 

LOWER WING 

FWD 



AFT 

6.68 

UPPER FUSELAGE 

FWD 

.55 


AFT 

.51 

UPPER WING 

FWD 

1.34 


AFT 

.61 

SIDE FUSELAGE 

FWD 

.64 


CORRUGATED SURFACE PERCENT CHANGE 



6.68 

3. 

NOT 

EVALUATED 



5.94 

-1. 

NOT 

EVALUATED 



7.52 

13. 


.56 

2. 


.43 

-16. 


2.29 

71. 


.87 

43. 


.80 

25. 


TABLE 7-V COMPARISON OF SMOOTH BODY AND CORRUGATED PANEL INSULATION REQUIREMENTS 



insulation materials) at infrared wavelengths. Although this mode of heat transfer 
is presumably accounted for in measurements of thermal conductivity, the work of 
Hughes (Ref. 27) showed that significant errors can occur when consideration of this 
phenomenon is omitted. 2) Two-dimensional conduction effects in the insulation layer 
on the upper and side surfaces of the vehicle. This would be caused by the small 
insulation thickness relative to the corrugation height and the large temperature 
gradients along the surface of the corrugation. 3) The previously mentioned 
variation of the correlated average heating with the changing flow conditions of 
the trajectory. 

Although these refinements of the insulation sizing problem should surely be 
considered in any further study, the present state of development of the heating 
correlations and definition of the vehicle inviscid and boundary layer flow con- 
ditions does not warrant their inclusion at this time. Further, for the present 
purposes of gauging the approximate insulation penalty associated with corrugated 
surfaces, omission of these factors in not of great importance. 

Application of the Wind Tunnel Correlations at Flight Conditions 

Two difficulties were encountered with regard to evaluating the full form 
MRA correlations (Equations 1, 3, 5, 7) at the conditions predicted by the inviscid 
flow analysis for the upper and side surfaces of the vehicle. The first was a 
failure of the correlation for the flow reattachment angle to predict values with- 
in an acceptable range, i.e., 0.0< sina < sina max . The correlation typically 
resulted in values of sina greater than sina mflX and often predicted sina greater 
than 1.0. This was initially overcome by arbitrarily limiting the value to sina m ; 
however, even with this limit imposed, the correlation for peak corrugation heating 
also consistently resulted in unreasonably high values (h m /h.. on the order of 

iTiaX T p 

10-20). It is felt that both of these difficulties result from the attempt to 
apply the correlations to flow conditions far removed from those included in the 
experimental data, and the ensuing possibility that some of the parameters included 
in the full form correlations which are of small statistical importance (i.e., 
poorly correlated at experimental conditions) result in large errors in the pre- 
dictions made at flight conditions. 

Use of the simplified forms of the MRA correlations (Equations 2 and 4) appear 
to avoid these problems and were for this reason finally used to obtain the 
corrugation heating for the leeward and side points of the vehicle. Since the 
full form correlations were used to obtain the corrugation heating at windward 
points, checks were made at the lower body forward and lower wing forward locations 
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using both forms of the correlations. Peak heating values calculated using the 
full form correlations were found to deviate only +5 to -17 % from those predicted 
by the simplified form. 

While the MRA correlations based on all available data (Equations 3 and 4) 
were used to determine the corrugation peak heating and thereby the peak temper- 
atures for skin material selection, the correlation based only upon thick boundary 
layer data (Equation 2) and the associated value of average heating was also 
evaluated for each leeward case investigated. The peak heating values based on 
the correlation of thick boundary layer data showed an average absolute deviation 
of 12% from those predicted by the correlation of all data. An exception was the 
value predicted for the upper wing forward location which was 47% lower than the 
value predicted by the correlation of all data. 

Corrugated Panel TPS Applicability 

It is apparent from the present analysis that there exist large areas on the 
SSTO for which use of the corrugated panel TPS is not appropriate. These include, 
as seen in Table 7- III, the lower wing at 20% chord for which the predicted peak 
corrugation temperature exceeds the 1644 K (2500°F) upper use limit of FS-85 
Columbium and may be presumed to include any other windward surface with relatively 
thin local boundary layer, i.e., leading edges or near the vehicle nose. Even the 
somewhat less severe case of the lower wing at 60% chord results in an increase in 
predicted peak temperature that spans the usable range of TD-Ni20C r . It is further 
apparent that even in windward areas for which the boundary layer is thick such as 
the lower fuselage at either the forward or aft points, peak corrugation temperatures 
dictate a change of material from that required for the smooth body case (typically 
from L-605 to TD-Ni20C r ) for all cases investigated. It may be generally concluded 
that although the use of corrugated panels on windward surfaces is not precluded by 
temperatures that exceed the use limits of currently available materials, the re- 
sulting peak temperature increments are typically of sufficient magnitude to re- 
quire material changes. 

The leeward and side fuselage points, although subjected to corrugation peak 
heating multipliers far greater than those of windward side locations (primarily a 
consequence of the high local flow angles present) experience significantly lower 
peak temperatures as a result of the much lower leeside smooth body heat transfer 
and recovery temperature. Material requirements are similarly reduced. With the 
exception of the upper wing at 20% chord, the predicted peak corrugation temperatures 
at all leeward and side locations investigated fall within the use limits of 6AL-4V 
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Titanium. The higher peak temperatures present at the forward wing location are, 
again, a result of the relatively thin boundary layer and consequent high smooth 
body and corrugation peak heating rates. It is concluded that the corrugated 
panel TPS should find wide application for leeward and side fuselage locations 
and is especially attractive since peak temperatures allow the use of the low 
density Titanium skin. 

As was previously pointed out in Table 7-V, departures of insulation thickness 
required for corrugated panels at the lower fuselage locations from those for the 
smooth body case are quite small (3% and -1%) and are, in fact, less than the limits 
of resolution of the present analysis. Application of the corrugated panels at 
upper fuselage locations results in similar small changes in absolute insulation 
thickness. Only two of the locations investigated, the lower wing at 60% chord and 
upper wing at 20% chord showed a significant insulation penalty for corrugated 
panels. The increases for these two points, 13% and 71% respectively, both 

o 

correspond to roughly 0.9 cm of added insulation thickness (56.1 kg/m Microquartz) 
or only about 45 kg per wing for either location. An even more moderate increase 
(0.16 cm) is predicted for the side fuselage location. The overall insulation 
penalties are by these estimates insignificant at fuselage locations aft of 50% 
of the vehicle length and of only small weight impact even in areas such as the 
lower wing for which peak corrugation temperatures exceed the 1478 K use limit of 
TD-Ni20C r . 
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8.0 AERODYNAMIC FLIGHT ASSESSMENT 

The change in entry range of the SSTO due to the drag associated with wavy 
wall panels was determined and the substantiating analysis is described in the 
following paragraphs. 

The nominal altitude, velocity, angle of attack and bank angle history for 
entry from 185 Km (100 nmi ) orbit is shown in Figure 7-2. For the first 3000 seconds 
of the entry, the angle of attack exceeds 30 degrees and the drag-due-to-1 i ft will 
dominate the drag. Changes in zero-angle-of-attack drag, Cd 0 , will have a small 
effect on total C^. This can be seen from Figure 8-1 which shows drag coefficients 
measured on a M scale model of the SSTO in the Langley Research wind tunnels 
(Ref. 28). At 30 degrees angle of attack the C^ is approximately 0.38; Cq q is 0.04 
at low angle of attack. A 10 percent change in C Dq results in only a 1 percent 
change in C Q . After 5000 seconds into the entry, the angle of attack is reduced 
to 10 deg. and the bank angle to 0 degrees. At these flight conditions, the C Dq 
is important. The wind tunnel data indicate C Do is 70 percent of the total C Q 
at Mach 4.63. 

The SSTO lift and drag coefficients were computed as a function of time from 
the trajectory data of Figure 7-2 (from Ref. (28)) by assuming the two-dimensional 
flat earth point mass equations of motion applied. The resulting total C Q and C^ 
are shown in Figure 8-2 as a function of altitude. The Mach regime of the wavy 
wall experiments corresponds to the altitude range from 27 to 39 Km (90 to 127 kft). 
Below this altitude range, the panel drag characteristics must be extrapolated. 

To estimate the effect of panel on range, the worst conditions for panel 
C n were selected. These conditions are: 

U /- 

Reynolds number = 10 x 1 0 /Meter 

Corrugation angle = 90 deg (perpendicular to freestream) 

Boundary Layer thickness, <$ = 2.5 cm 

The experimental panel drag coefficient data shown in Figure 8-3 give a C D of 
0.0107 at the Mach 2.4 worst case condition. To estimate the worst case Mach 
4.5 C^, the data in Figure 8-4 for the thick boundary layer were used. The ratio 
for a thin and thick boundary layer at Mach 2.4 is 0.0107/0.0085 = 1.26. The 
thick boundary layer at Mach 4.5, 0.0036, was multiplied by 1.26 to give a 
of 0.0045 for the thin boundary layer at Mach 4.5. An inverse dependency of 
with Mach number was assumed to provide an estimated at Mach 10 of 0.0020. 

This 1/Mach dependence is slightly weaker than the dependency indicated in Figure 
9-4 which is (1/Mach)l .37. 
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These individual panel C^'s were related to the total SSTO drag coefficient 
using the following expression 


aC D C d X S 


panel 


ref 


total 

5 panel 


where S 


2 

panel = ex P enmental panel area for (8.01 cm ) 


5 ref 


total 


= SSTO reference area 881 (9484 ft^) 

= total panel area on SSTO (parametrically varied) 


By varying the maximum change in SSTO drag coefficient, aC d , can be determined 

Figure 8-5 summarizes these results for Mach 2.4, 4.5 and 10.0. The aC d varies 

linearly with total panel area. For the entire SSTO paneled, maximum aC^'s of 0.034, 

0.014, and 0.006 can result at Mach 2.4, 4.5 and *10.0 respectively. A more practical 

2 2 

panelling area would correspond to the leeside body area of 703 m (7571 ft ) which 
results in a aC^ of 0.0085, 0.0035, and 0.0016 for the same Mach numbers. 

In Figure 8-6 these data are compared with the computed from the trajectory 
data. For Mach numbers greater than 4.5 the aC^ is small compared to the total 
drag coefficient and probably can be neglected. However, for Mach numbers less 
than 4.5, the aC^ increases substantially. For a fully panelled SSTO the aC^ 
approximates C Q at Mach 2.0 if a 1/Mach dependency occurs. For the leeside body 
panelling only, the relative magnitudes are reduced but become large at Mach 1.0. 

These results were related to trajectory range change by perturbing the C D of 
the basic trajectory by a aC^ and computing the resultant change in range from 
entry to a given altitude. The results are shown in Figure 8-7. From entry to 
30Km (100 kft) altitude the range change is a maximum of 24Km (13 nmi ) for a worst 
case, fully panelled SSTO. Below this altitude the Mach number is reduced and the 
aC q is a larger contribution to C Q . As a result the range change from entry to 
9Km (30 kft) can be appreciable. For the fully panelled SSTO, a 83Km (45 nmi) 
range loss can occur. For leeside body only, the range loss is 24Km (13 nmi). 
Insufficient panel data and trajectory information exist below 9Km (30 kft) 
to assess the full impact at these conditions. However, increased C D due to wavy 
wall effects could substantially reduce range below 9Km (30 kft) altitude. (Note: 

The aCq for Mach numbers below 2.4 was held constant for the Figure 8-7 results and 
could be appreciable higher if the 1/Mach dependency applies). 

These data were cross-plotted as a function of panel area in Figure 8-8. 

The variation of range with panel area is readily apparent from the cross 
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plot. The conclusions from these results are as follows. Wavy wall drag in- 
crease has little impact on the SSTO entry trajectory above 30Km (100 kft) 
altitude even in the worst case. Below this altitude, panels oriented perpendicular 
to the freestream can have a significant effect on range. Below 9Km (30 kft), in- 
sufficient data exist to assess the impact on range. However, the trend in the 
panel Cj variation with Mach number appears to be an inverse relationship which 
could result in large range penalties at the final approach speeds. More data 
are required to quantify this low speed effect. 
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9.0 STRENGTH AND FLUTTER ASSESSMENT 

For the present corrugated panel flight evaluation, a single corrugation 
shape (Figure 7-5) based upon longitudinal bending stiffness requirements and 
previous experience with the Gemini TPS was assumed throughout the analysis. 

No attempt was made to optimize the panel shape or thickness as a function of the 
final calculated temperatures. The strength assessment yielded a panel length 
(distance between supports) ranging from 20.3 to 50.8 cm and temperature use 
limits for the corrugation materials considered as follows: 

Columbium (FS-85) 1644 K (2500°F) 

TD Ni - 20 1478 K (2200°F) 

L-605 1256 K (1800°F) 

Rene' 41 1144 K (1600°F) 

Titanium (6AL-4V) 811 K (1000°F) 

To determine the suitability of such a panel shape and these material 
selections for application in the SSTO dynamic environment, a flutter analysis was 
performed for the flow normal to the corrugations using the design approach given 
by Ref. 29. 

Results of this analysis indicated that flutter does not occur for any of the 
cases studied for panel lengths up to the 50.8 cm derived from the strength con- 
siderations. Since the flow was assumed to be normal to the waves, flutter should 
not be a problem for other local flow angles for the high Mach number cases con- 
sidered in this study. 
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10.0 SUMMARY OF RESULTS 

A wide range of wind tunnel data for aerodynamic heating and pressure 
distributions on corrugation roughened surfaces in thick supersonic and hyper- 
sonic turbulent boundary layers has been analyzed and correlated. The effects of 
the corrugations in terms of increased aerodynamic heating and drag and the 
initiation of flutter were also considered for an advanced space transportation 
(a single-stage-to-orbit) vehicle. The data analysis and flight assessment yielded 
the following results: 

1. For almost all experimental cases studied, the flow separated on the 
expansion surface and reattached on the compression surface. 

2. Peak heat transfer increased as the angle between the corrugations 
and the free stream flow increased, except for one experiment at the 
lowest Reynolds number (Re^/M = 1.3 x 10 ). 

3. The peak heat transfer rates showed a direct proportionality to Reynolds 
number except for the deepest bead tested at the highest Reynolds 
numbers in which cases the peak heating was inversely proportional to 
Reynolds number. 

4. Heat transfer and pressure were fairly insensitive to wave amplitude unless 
the corrugation protruded far out into the boundary layer. Severe heating 
increases were observed for e/&*~ 0.1 in hypersonic flow. 

5. Changes in wavelength and surface wave radius by a factor of two produced 
a small effect on heating. 

6. For some flow conditions, the average wave heating was less than the 
corresponding flat plate value. The functional relation between average 
heating and corrugation angle was not consistent with maximum heating 
trends in that the average heating decreased with increasing flow angle 
for many conditions. 

7. The measured pressure drag coefficients on corrugated panels were extremely 

-7 -3 

low for hypersonic flow, ranging from 8.3 x 10 to 2.16 x 10 . 

8. Peak pressure measured on the compression surface could not be explained 
on the basis of the local internal flow passing through an imbedded shock 
wave. For the flow normal to the corrugations, the maximum pressure 
correlated in terms of the total pressure and Mach number of the stream- 
line at the wave height. 
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9. Shifts in the peak pressure location were detected. Based on previously 
obtained data for flow over notches, the changes in the peak pressure 
location could be explained by a separated flow region dependent on the 
Mach number at the wave height. 

10. Peak heating values for all available data for thin and thick turbulent 
boundary layers were correlated in terms of bulk boundary layer, internal 
boundary layer, and geometric parameters. 

11. Average heating data for thick turbulent boundary layers were correlated 
in terms of the corresponding maximum heatinq values, the local flow con- 
ditions, and the geometric parameters. 

12. The geometry of the separated flow region was correlated in terms of the 
local surface angle at the point of maximum heating. This new formulation 
improved the capability of the heating correlations. 

13. For the thermal flight assessment, severe heating penalties were found to 
exist over large areas of the windward surfaces with thin boundary layers, 
leading edges, and locations near the nose. Temperature penalties ranged 
from 100K to 500K, and very large structural advantages of the corrugations 
would be required to overide these heating penalties. 

14. The use of corrugated panels on windward surfaces with thick turbulent 
boundary layers could dictate a change in material from that required for 
the smooth body case. 

15. Computed peak corrugation temperatures were within allowable limits of the 
equivalent smooth wall material for leeward and side fuselage locations. 
Hence, corrugated panel TPS have application for these areas. 

16. Insulation weight penalties caused by the use of corrugated panels were 
insignificant for all vehicle locations where the panels were found to 
have application. 

17. The change in entry range of the SSTO due to increased drag associated 
with corrugated panels was small for altitudes above 30Km (100K ft). For 
a fully panelled SSTO, a maximum reduction in range of almost 10 percent 
was predicted from entry to 9Km (30K ft) altitude. 

18. Additional corrugation drag data are needed below Mach 2 to assess the 
impact of corrugations on range at final approach speeds. 

19. Panel flutter was not predicted for the SSTO flight trajectory in the 
supersonic and hypersonic flow regimes. 
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12.0 LIST OF SYMBOLS 


'DO 


"f 

CFHT 


h 

k e 

L 

M 

P 

I 

P t 

q 

Re 

co 

Re 

e x 

Re 

e 

Re ,Rg 

. e 

r 

S 

SSTO 

T 

UPWT 

U 

X 


speed of sound 
SSTO drag coefficient 


2 x DRAG 
pV ! S REF 


C Q at zero angle of attack 
Panel drag coefficient 
skin friction coefficient 
Continuous Flow Hypersonic Tunnel 
hydraulic diameter 

heat transfer coefficient, q/(T - T ); also enthalpy 

equivalent sand-grain roughness height 

wavelength 

Mach number 

pressure 

Pitot pressure 

heat transfer rate 

p U 

CO CO 


free-stream Reynolds number, 

°° y x 

Reynolds number based on surface length, p °° °° 


Reynolds number based on wave amplitude, 


^00 

p U e 


u 


p U e 

Reynolds number evaluated at wave amplitude, £ - - £ — 


surface wave radius 

surface distance measured from top of wave; reference area 

Si ngl e-Stage-to-Orbi t-Vehi cl e 

temperature 

Unitary Plan Wind Tunnel 
velocity 

axial distance measured from top of wave; smooth wall surface distance 
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sep 

Y 

Z 

a 


y 

e 


6 * 


6 

e 

u 


p 


distance to separation point 

spanwise distance along tunnel sidewall; also wave vertical coordinate 

skin thickness; also distance through boundary layer 

local wave surface angle 

ratio of specific heats 

maximum wave amplitude from wave midline 

laminar sublayer thickness 

displacement thickness 

boundary layer thickness 

momentum thickness 

molecular viscosity 

angle of corrugations relative to free-stream flow direction 
(See Figure 4-1) 

density 
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I 


SUBSCRIPTS 


av 


aw 


EQ 


( ) 

( ) 

( ) 

( ) 

< >1 
( ) 

( ) r 

l >t 

^ ^ TOTAL 

( ), 


FP 


MAX 


'w 


( ), 


( )« 
( )• 


average 
adiabatic wall 
equil ibrium 
flat plate 
local 
maximum 

reattachment point taken as point of maximum pressure 
stagnation conditions 
total distance 

wall conditions (same as flat plate conditions) 

evaluated at maximum wave amplitude 
free-stream conditions at edge of boundary layer 
in front of shock 
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APPENDIX A 


WAVE SHAPE FOLD OUT 



FIGURE A-l CORRUGATION SHAPES INVESTIGATED 
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